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THE  EFFECT  OF  THERaOPECRblS  OK  THE  PROCEED  OF  PRECIPITATION 
OF  AdH  PARTICLE-*  CN  COGnED  BLALEj  OF  GAo  TURBINES 


S.  Sh.  Rozenberg,  Engineer 


Certain  ideas  are  presented  concerning  the  important 
influence  exerted  by  the  thermophoreais  effect  on  the 
process  of  precipitation  of  particles  on  cooled  blades 
of  gas  turbines.  Experiments  on  a  cooled  and  an  uncooled 
plate  (a  blade  model)  showed  that  in  comparison  with 
other  factors  this  phenomenon  exerts  the  most  important 
influence  on  the  process  of  precipitation  of  particles. 


Modern  gas-turbine  construction  is  characterized  by  the  tendency  to  increase 
the  gas  temperature  in  front  of  the  turbine  while  cooling  the  most  heavily  stressed 
components  of  the  machine  and  by  the  tendency  to  use  inexpensive  types  of  fuel. 
However,  owing  to  the  presence  of  ash  in  the  combustion  products  of  heavy  liquid 
and  solid  fuels,  the  use  of  these  fuels  causes  additional  complications  resulting 
from  erosion,  corrosion,  and  precipitation  of  ash  in  the  flow-through  sections  of 
gas  turbines.  Practice  shows  that  these  difficulties  increase  in  high-temperature 
turbines.  Thus  the  study  of  the  process  of  precipitation  on  cooled  surfaces, 
including  blades,  is  of  indubitable  interest  in  the  field  of  gas-turbine 
construction. 

The  precipitation  of  ash  on  turbine  blades  is  a  complex  physicochemical 
process,  in  which  forces  of  different  origin  may  play  a  role,  depending  on  the 
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conditions  of  precipitation  of  the  solid  particles*  These  forces  may  be  mechanical 
hydrodynamic,  chemical,  molecular,  thermal,  or  electrostatic.  One  of  the  causes 
of  precipitation  of  ash  particles  and  their  retention  on  a  surface  is  the  action, 
of  thermophoresis  forces*  is  is  known,  the  thermophoresis  effect  occurs  when 
aerosol  particles  suspended  in  a  temperature  field  move  into  a  region  of  lower 
tempera t urea,  i.e*,  towards  a  cold  body  [ll* 

The  thermophoresis  force  acting  on  a  spherical  aerosol  particle  in  the 
absence  of  convection  currents  without  taking  into  account  radiative  heat  transfer 
from  the  particle  can  be  calculated  from  the  formula  given  in  Ref.  [ljt 


F,  =.-  -  9n 


ha 

2  Xa  X/ 


jSr 
QT ‘ 


dT 

dx 


U) 


where  F ^  is  the  tbermophoresis  force; 

r  is  the  radius  of  the  particle; 

M  and  f  are  the  viscosity  and  density  of  the  gas; 

T  is  the  absolute  temperature  of  the  gas; 
and  are  the  thermal  conductivities  of  the  ga3  and  the  particle; 
dT/dx  is  the  temperature  gradient  in  the  gas  along  the  x-axis. 

In  the  case  of  uniform  temperature  distribution  in  an  incident  flow  passing 
around  a  surface  the  temperature  of  which  is  different  from  that  of  the  flow  the 
maximum  temperature  gradient  will  occur  in  the  thermal  boundary  layer  [2].  Since 
the  thermophoresis  force  is  directly  proportional  to  the  temperature  gradient,  it 
follows  that  the  thermophoresis  effect  will  be  most  strongly  pronounced  in  the 
boundary  layer,  where  it  can  have  a  decisive  effect  on  the  nature  of  the  motion 
of  the  particle. 

In  particular,  if  the  temperature  of  the  wall  is  lower  than  that  of  the  flow, 
e*g*,  in  the  case  of  cooled  turbine  blades,  the  thermophoresis  force  will  cause 

the  trajectory  of  motion  of  the  particle  to  be  deflected  towards  the  wall.  This 

may  considerably  increase  the  number  of  particles  colliding  with  the  wall  or 
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precipitating  on  it,  in  comparison  with  the  case  of  an  uncooled  surface,  where 
the  temperature  of  the  flow  and  the  wall  is  the  same  and  the  thermophoresia  effect 
is  absent. 

Since  the  motion  of  a  particle  about  a  turbine  blade  depends  on  various  forces, 
in  order  to  determine  the  effect  of  thermophoresis,  it  is  necessary  to  schematize 
to  phenomena .  Let  us  consider  the  motion  of  particles  in  the  laminar  boundary 
layer  near  a  flat  wall.  This  will  enable  us  to  ascertain  the  basic  tendencies  of 
the  given  phenomenon. 

In  setting  up  the  equations  of  motion  of  a  particle  in  the  laminar  boundary 
layer,  let  us  assume  (fig.  l)  the  following! 

1.  The  thermophoresis  force  is  calculated  from  formula  (l),  the  value  of  T 
in  this  formula  being  taken  equal  to  the  average  temperature  of  the  boundary  layer 

T  —  r_7'o  +  7'*f 

2 — •  (1‘) 

2.  All  the  physical  constants  are  assumed  to  be  constant  for  the  entire 
thickness  of  the  boundary  layer  and  are  determined  for  the  temperature 

/  -  ‘*+‘o  Cl") 

ST  2  • 

3.  For  the  temperature  distribution  in  the  boundary  layer  we  take  the 
dependence 


f  -  9.(2%  -  %)’• 


(2) 


Here 


(3) 

(3‘) 
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where  T0»  T»  and  Zy  are  the  absolute  temperatures  of  the  flow  outside  the  boundary 

layer,  within  the  boundary  layer,  and  on  the  wail,  respectively! 


’I  T  is  the  dimensionless  ordinate  and 


(4) 


Ug.  1*  Diagram  for  determining  the 
trajectory  in  the  boundary  layer* 


The  thickness  of  the  thermal  boundary  layer  b^,  as  is  known,  is  related  to 
the  thickness  of  the  hydrodynamic  boundary  layer  by  the  dependence  [31 


«T=*Pr  3.  (5) 

Taking  Pr^Q«7  for  air,  we  have 

«T  (5.) 

The  thickness  of  the  hydrodynamic  boundary  layer  is  calculated  from  the 


formula  [4l* 


6  =  5,83  ]/  =  5,83 

f  C0 


(6) 


where  *  is  the  kinematic-viscosity  factor  of  the  gas. 

4*  Tor  the  distribution  of  the  velocity  component  c  we  take  the 


approximation 


Here 


=  Co  (2n  —  ii1). 


tl  =  T: 


(7) 


(6) 
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«q  is  the  flow  Telocity  outside  the  boundary  layer) 

*  is  the  dimensionless  ordinate  of  the  hydrodynamic  boundary  layer. 

5*  Tor  the  distribution  within  the  boundary  layer  of  the  velocity  component 
c  perpendicular  to  the  wall  we  take  an  approximate  dependence  obtained  by  the 

y 

author  by  solving  the  continuity  equation  with  the  variation  of  the  gas  temperature, 
and  consequently  the  variation  of  the  density  in  the  boundary  layer,  taken  into 
account* 

This  dependence  has  the  form 


(9) 


or 


5,83v  ,  ,  . 

c„  =  — r—  (an*  — 


Let  us  introduce  the  dimensionless  temperature  parameter  k 

T  T* 

*  -  Is.  = . 


(9‘) 


(10) 


”u  Ta  rCT 

The  coefficients  £  and  are  functions  of  this  parameter. 

Turning  to  the  derivation  of  the  equations  of  motion  of  a  particle  in  a 
nonisothermal  boundary  layer,  let  us  consider  the  action  of  the  forces  of  inertia, 
the  resistance  of  the  medium,  and  thermophoresis. 

Preliminary  calculations  show  that  the  velocities  of  the  particles  relative 
to  the  flow  are  such  (small  Re  numbers)  that  we  can  use  Stoke's  law  [5]  to  express 
the  resistance  of  the  flow  to  the  motion  of  the  particle. 

Then,  assuming  the  motion  to  be  planar,  we  have 

tlj  d*y  _  _  . 

m  rfT  ■-=  m  Hi*  =  F>v  ~  <1' — cvy' 


dt 

du 


d*x 


m~df  =  m  Hi*  =  F,x  ~  &nVLr  (u  ~  Ct)' 


(11) 


where  x  and  ji  are  the  projections  of  the  velocity  of  the  particle  on  the  y-axis 
and  x-axis,  respectively; 

r  and  are  the  projections  of  the  thermophoresis  force  determined  from  formula 
ty  tx 


HD-TP-62-1025/1+2H 


■5 


(l)  on  the  y-axis  and  the  x-axis { 

o 

A:  is  the  mass  of  the  particle. 

From  equations  (2)  and  (4)  we  have 


—  =  —  M1  -  ’It)// 
dx  Xbt 


From  formulas  (11')  and  (11")  it  can  be  seen  that 


<96 


db 


dx  ^  dy 


Hence,  taking  into  account  formula  (l),  we  obtain 

Fu  «  F,y. 


(ll*) 

Ui") 


(12) 


(12*) 


This  enables  us  to  neglect  the  longitudinal  component  of  the  thermophoresis  force 
in  equations  (ll). 

Next,  let  us  consider  the  two  cases,  where  the  particles  l)  do  not  hare 
considerable  velocities  and  2)  have  a  considerable  velocity  relative  to  the  flov. 
Of  practical  interest  is  the  case  where  the  particle  velocity  is  directed  towards 
the  wall. 

If  the  particle  does  not  have  a  considerable  velocity  relative  to  the  flow 
and,  as  a  consequence  of  this,  does  not  undergo  important  accelerations  or 
retardations  resulting  from  the  resistance  of  the  medium,  we  can  neglect  the 
inertial  terms,  and  then  equations  (ll)  with  (l)  and  (12')  taken  into  account  will 
assume  the  following  formi 


-9.1 


21,  +  ki 


6g  (2  —  2r|r)  4 - 

or 


—  6  (v  —  cy)  =  0; 


u  —cx  =  0.  j 

After  transformations  we  obtain  the  system  of  equations 


(13) 


«  =  c,  (2r|  —  n*). 


CM) 
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Let  us  introduce  the  dimensionless  coordinates  x^  and 


*  =  -f 

0» 


(15) 

(15') 


where  xQ  is  the  coordinate  of  the  entry  of  the  particle  into  the  hydrodynamic 
boundary  layer; 

t  is  the  thickness  of  the  hydrodynamic  boundary  layer  at  this  spot. 
Taking  (6)  and  (8)  into  account,  we  have 

yx  *i  Vx 

Then,  letting 

P 


(16) 


3  Qo  pr^  **  _ 
5,83*  T  2  ka  +  X, 


(17) 


and  bearing  in  mind  that  y  and  =  u,  we  obtain  after  a  number  of 
transformations  and  with  (lo)  taken  into  account 


dx. 


L-(a»l*  — />!)*)  —  ?  (pr  3  -  r,) 


5.83 


2r|  —  tj* 


2L 
2  • 


(18) 


This  differential  equation  is  solved  by  assigning  values  to  the  coefficients 
£  and  which,  as  was  noted  above,  are  functions  of  the  temperature  parameter 
In  particular,  solving  this  equation  for  Pr  *  0.7  and  k  =  2  with 

b  =  1,644.  |  (18') 


and  taking  into  account  the  initial  conditions 

•*i  -  i ;  *i =  i •  (18”) 

we  obtain  a  formula  relating  the  coordinates  of  a  particle  moving  in  a  boundary 

layer  in  the  presence  of  thermophoresis  forces 

-dWH-MW  ■-[2,7sSn._  + 


+  (0,409  +  1 ,6P)  In  I  ■24y55~6~T)  1  +  ^J  = 

Kt-0-13 


1,37  V?  —  0,13?* 


arCban 

1,37 


p 

T)  -  0,5p 


•0.5P  1 

-o,i3p*  y 


(19) 
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ngur*  2a  shows,  in  the  coordinates  xL  and  the  trajectories  of  a  particle 
for  different  values  of  the  parame ter  0  with  k  =  2.  From  the  diagram  it  can  be 
seen  that  the  trajectory  of  motion  of  the  particle  is  strongly  dependent  on  the 
value  of  0.  The  greater  the  value  of  the  greater  the  thermophoresis  force  end 
the  more  rapidly  the  particle  reaches  the  wall.  Thus,  under  the  action  of  a 
thermophoresis  force,  even  a  particle  that  does  not  have  a  velocity  component  in 
the  direction  of  the  wall  at  the  initial  moment  will  atill  strike  the  wall, 
assuming,  of  course,  that  the  wall  is  sufficiently  long. 

If  the  particle,  upon  entering  the  boundary  layer,  has  a  considerable 
velocity  relative  to  the  laminar  boundary  layer,  resulting,  for  example,  from 
pulsations  and  vorticity  in  the  incident  flow  or  from  the  action  of  centrifugal 
forces,  we  can  no  longer  neglect  the  inertial  terms. 

Of  practical  interest  is  the  case  where  the  velocity  of  the  particle  is 
directed  perpendicularly  to  the  wall,  for  in  this  case  its  velocity  may  be 
comparable  to  the  velocity  component  of  the  flow  in  the  direction  of  the  wall  c 

y 

and  to  the  velocity  arising  as  a  result  of  thermophoresis. 

In  deriving  the  equations  of  motion,  let  us  neglect  the  inertial  term  along 
the  x-axis  and  take  it  into  account  along  the  y-axis. 

Then  equations  (li)  will  assume  the  form 

dv  d*u 

OT  flfT  =  m~dt*~F,v-  6n»‘r  -  cv)< 

u  —  cx  =  0.  .  . 

(20) 

Switching  over  to  dimensionless  coordinates,  after  a  number  of  transformations 
we  obtain  a  syateju  of  equations  describing  the  motion  of  a  particle  in  a  laminar 
nonisothermal  boundary  layer  for  Pr  ^  1  in  the  presence  of  an  initial  velocity  of 
the  particle  in  the  direction  of  the  wall 


where 


r 


wi 


(21) 

V 

' 

t 

(22) 

X  ’ 

c\j 

.  -!i£i 
vs  ’ 

(M 

Here  is  the  specific  gravity  of  the  particle. 


Fig.  2.  The  path  of  the  particle  iu  the  direction 
of  the  wall  in  the  laminar  boundary  layers  a)  with¬ 
out  initial  velocity;  b)  with  initial  velocity!  l) 

0  »  0.02;  2)  0  *  0,01;  j)  with  thermophoresia  forces 
taken  into  account;  4)  without  taking  thermophoresia 
forces  ihto  account. 


The  quantity  is  called  the  "relaxation  time"  and  is  widely  used  in  aerosol 
mechanics  [lj. 

The  system  of  equations  (2l)  cannot  be  solved  in  the  general  form  and  must  be 
solved  numerically  in  each  particular  case.  In  order  to  estimate  the  nature  of  the 
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motion  of  a  particle  with  a  certain  initial  velocity  in  the  direction  of  the  wall, 

we  uaed  Adam1 a  method  [6]  to  obtain  a  numerical  solution  to  the  system  of 

equations  (2l)  for  the  following  initial  conditional 

Tq  *  973°K  3  *  0.02 

1c  s  2  He  ■  511 

0 

p  *  1  atm  aba  Tj_  3  *0.06 

cQ  *  20  m/aee  Yj  *  2200  *8/*°? 

r  »  5  .  10*6  m  (d  »  10  m)  80  3  5.16  .  10-4  m 

The  particle  trajectory  calculated  from  the  system  of  equations  (21),  which 

take  into  account  the  action  of  thermophoresis  forces,  is  shown  by  curve  )  in  Fig. 
2b.  Curve  4  shows  the  trajectory  of  the  same  particle,  but  without  taking  into 
account  the  action  of  thermophoresis  forces.  As  can  be  seen  from  the  diagram, 
in  the  first  case  the  particle  strikes  the  wall,  while  in  the  second  case  the 
particle,  having  'pierced*  as  a  result  of  its  kinetic  energy  part  of  the  thickness 
of  the  Doundary  layer,  after  its  velocity  in  the  direction  of  the  wall  is  attenuated 
will  follow  the  line  of  the  gas  flow  and  gradually  veer  away  from  the  wall.  This 
clearly  shows  that  in  the  case  of  a  particle  moving  near  a  cooled  wall  and  possess¬ 
ing  an  initial  velocity  in  the  direction  of  this  wall  the  phenomenon  of 
thermophoresis  also  has  a  decisive  effect  oa  nature  of  the  motion  of  the 
particle. 

For  a  given  thermal  conductivity  of  the  flow  the  thermal  conductivity  ^ 

of  the  particle  itself  has  a  strong  effect  on  the  nature  of  the  motion  of  the 
particle.  Obviously,  the  thermal  conductivity  of  the  particle  is  determined  by 
its  chemical  «pd  mineralogical  composition,  and  also  by  its  structure,  by  the 
presence  of  pores,  cavities,  etc.  Moreover,  the  velocity  of  motion  of  a  particle 
under  the  action  of  tbermophoresis  is  affected  by  its  shape,  which  is  also  dependent 
on  its  mineralogical  composition  and  structure.  Thus,  it  is  clesr  that  particles 
of  different  chemical  and  mineralogical  composition  and  different  crystalline 
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structure  vill  have  different  trajectories  In  the  boundary  layer  near  the  cooled 
wall  and  will  move  and  strike  the  cooled  well  in  different  ways,  as  a  result  of 
which  particles  with  a  lower  thermal  conductivity  will  approach  the  wall  and 
precipitate  on  it  faster.  This  explains  the  well-known  fact  [7]  that  different 
particles  precipitate  in  different  ways  on  a  wall  that  is  colder  than  the  flow. 
It  is  known  that  in  this  case  the  particles  that  precipitate  on  the  wall  are 
primarily  particles  of  definite  composition,  which  may  be  due  to  their  differing 
behavior  in  a  temperature  field  under  the  action  of  thermophoresis. 

Analyzing  the  equations  of  motion  of  a  particle  in  a  laminar  nonisothermal 
boundary  layer,  it  is  obvious  that  the  motion  of  a  particle  with  an  initial 

velocity  will  be  determined  by  the  following  dimensionless  complexes i 

2 


n-  ■  R  —  _JL_  .  >'u 

«ex,  -  v  .  P  -  5  83,  T  2>.u  - 

J*.  ■  „  v»  ■  b  T» 

^ «  % 


PrJ 


Note  that  the  quantity  c^r  is  called  the  'inertial  path"  of  a  particle  with 

an  initial  velocity  cQ  [l]«  it  is  equal  to  the  path  traversed  by  a  particle 

with  an  initial  velocity  c_  in  a  medium  devoid  of  external  forces.  Thus  the 
rc  u 

complex  -y*  characterizes  the  'piercing  power'  of  the  particle  with  reference  to 
the  thickness  of  the  boundary  layer. 

Special  mention  should  be  made  of  the  complex  3  and  the  temperature 
parameter  k,  which  in  the  case  v^  =  0  completely  determine  the  motion  of  the 
particle  in  the  dimensionless  coordinates  x^,  y^.  In  this  case  they  are  the 
criteria  determining  the  nature  of  the  motion  of  the  particle.  The  temperature 
parameter  It  characterizes  the  intensity  of  the  cooling  of  the  wall,  while  3 
characterizes  the  ratio,  between  the  thermophoresis  forces  and  the  hydrodynamic 
forces. 

for  the  purpose  of  making  a  qualitative  experimental  check  of  the  above- 
mentioned  assumptions  concerning  the  important  influence  of  the  thermophoresis 
effect  on  the  process  of  precipitation  of  particles  on  a  cooled  wall  we  carried 
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out  special  experiment#  with  plates,  on  the  surfaces  of  which  different 
temperatures  were  maintained  and  which  were  placed  in  a  hot  gas  flow  (t  *  6lO^C) 
containing  fine  particles  of  ash. 

Three  pairs  of  plates  were  set  up  in  the  flowt  l)  cooled  with  a  wall 
temperature  of  2)  not  cooled  with  a  wall  temperature  of  607^|  3)  heated 

with  a  wall  temperature  of  680^. 

Experiments  with  a  cooled  and  an  uncooled  plate  showed  that  the  cooled 
plate  wears  out  many  times  faster  than  the  uncooled  plate.  As  a  result  of  an 
experiment  with  a  heated  and  an  unheated  plate  it  was  ascertained  that  precipitatioi 
was  almost  entirely  absent  on  the  heated  plate.  £mall  depositions  occurred  only 
at  the  sites  of  roughnesses  and  aerodynamic  shadow.  Obviously,  in  this  case  the 
conditions  of  precipitation  of  ash  particles  on  the  plate  surfaces  were 
identical  in  every  respect,  except  the  direction  of  the  heat  flux.  Thus  the  great 
difference  in  the  amount  of  precipitation  is  due  to  the  fact  that  in  the  ease  of 
the  heated  plate  the  thermophoresis  forces  were  directed  away  from  the  plate  and 
prevented  ash  particles  from  precipitating  on  the  surface.  The  small  depositions 
that  occurred  in  the  zone  of  aerodynamic  shadow  and  breakaway  may  be  due  to  the 
fact  that  particles  lanaing  in  the  vortices  of  the  flow  received  therefrom  a 
velocity  in  the  direction  of  the  wall  greater  than  the  velocity  received  from  the 
thermophoresis  force  counteracting  this  motion. 

Obviously,  in  the  experiment  with  the  cooled  plate  the  considerable 
depositions  on  it  are  the  consequence  of  the  thermophoresis  forces  acting  on  the 
particles.  These  forces  are  airectea  toward  the  wall  and  further  the  increase 
of  depositions.  Thus  the  experiments  carried  out  corroborate  the  theory  stated 
above  concerning  the  important  influence  of  the  thermophoresis  effect  on  the 
process  of  precipitation  of  aerosol  particles  on  cooled  or  heated  walls. 

On  the  basis  of  the  aforesaid  we  can  draw  the  following  conclusions* 

1.  In  the  formation  of  depositions  on  cooled  surfaces  the  action  of 
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thermophoresis  is  basic  in  comparison  with  other  factors  affecting  the  process 
of  precipitation.  This  is  of  particularly  great  importance  in  the  case  of  high- 
temperature  cooled  turbines  operating  on  combustion  products  containing  volatile 
ash.  In  the  case  of  gas  turbines  with  a  uniform  temperature  distribution  in  the 
incident  flow  the  thormophoresia  effect  will  manifest  itself  mainly  in  the 
boundary -layer  region,  where  the  maximum  temperature  gradient  eeeurs. 

2.  A  particle  impinging  on  a  nonisothermal  boundary  layer  near  the  cooled 
wall,  under  the  action  of  thermophoresis  forces,  will  necessarily  strike  the 
wall,  provided  that  the  wall  is  sufficiently  long. 

3.  The  nature  of  the  motion  of  a  particle  in  a  nonisothermal  boundary 

T  T 

layer  is  influenced  decisively  by  both  the  temperature  parameter  k  ~  .IS  ...» 

WQ  T0  "  TW 

and  by  the  ratio  between  the  thermal  conductivities  of  the  particle  and  the 
medium  A^/  A^. 

4.  By  choosing  the  system  of  coordinates  x.  ■  “  y  a  in  the  case  of 

1  x0  1  *  0 

absence  of  an  initial  velocity  of  the  particle  relative  to  the  flow,  the  nature  of 

the  motion  of  the  particle  in  a  laminar  boundary  layer  will  depend  only  on  the 

T 

temperature  parameter  k  *  ^  and  the  dimensionless  quantity 

*0  1 

r  _  3  .  j!o_  Pr3 

P_  5,832  T  2Aa+li  ‘ 

5.  In  the  case  of  greatest  practical  interest,  where  the  particle  has  an 
initial  velocity  in  tne  direction  of  the  wall,  the  motion  of  the  particle  in  the 
coordinates  x^,  y^  will  be  determined  by  the  complexes 


Re 


*0 


Vo  . 


P  = 


Pr3  Aq 


5,83s  T  2Xa  +  Xt  ' 
*  =  .!«-•  TCo 


t>i 


6.  Calculations  made  on  the  basis  of  the  above  theory  show  that  particles 
with  different  cnemical  and  mineralOfeical  compositions  and  different  thermal 
conductivities  and  shapes,  when  in  a  field  with  a  temperature  gradient,  must 
move  along  different  trajectories  with  the  particles  of  lower  thermal  conductivity 
striking  the  wall  sooner.  This  explains  the  well-known  phenomenon,  where 
deposition  of  different  ash  components  on  cooled  surfaces  occurs  dissimilarly 
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and  'selective  deposition*  takes  piece, 
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